QTQTN motif upstream of the furin-cleavage site plays a key role in SARS-CoV-2 infection and pathogenesis.

The furin cleavage site (FCS), an unusual feature in the SARS-CoV-2 spike protein, has been spotlighted as a factor key to facilitating infection and pathogenesis by increasing spike processing. Similarly, the QTQTN motif directly upstream of the FCS is also an unusual feature for group 2B coronaviruses (CoVs). The QTQTN deletion has consistently been observed in in vitro cultured virus stocks and some clinical isolates. To determine whether the QTQTN motif is critical to SARS-CoV-2 replication and pathogenesis, we generated a mutant deleting the QTQTN motif (ΔQTQTN). Here, we report that the QTQTN deletion attenuates viral replication in respiratory cells in vitro and attenuates disease in vivo. The deletion results in a shortened, more rigid peptide loop that contains the FCS and is less accessible to host proteases, such as TMPRSS2. Thus, the deletion reduced the efficiency of spike processing and attenuates SARS-CoV-2 infection. Importantly, the QTQTN motif also contains residues that are glycosylated, and disruption of its glycosylation also attenuates virus replication in a TMPRSS2-dependent manner. Together, our results reveal that three aspects of the S1/S2 cleavage site-the FCS, loop length, and glycosylation-are required for efficient SARS-CoV-2 replication and pathogenesis.

Evaluation of SARS-CoV-2 ORF7a Deletions from COVID-19-Positive Individuals and Its Impact on Virus Spread in Cell Culture.

The spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), causing the COVID-19 outbreak, posed a primary concern of public health worldwide. The most common changes in SARS-CoV-2 are single nucleotide substitutions, also reported insertions and deletions. This work investigates the presence of SARS-CoV-2 ORF7a deletions identified in COVID-19-positive individuals. Sequencing of SARS-CoV-2 complete genomes showed three different ORF7a size deletions (190-nt, 339-nt and 365-nt). Deletions were confirmed through Sanger sequencing. The ORF7a∆190 was detected in a group of five relatives with mild symptoms of COVID-19, and the ORF7a∆339 and ORF7a∆365 in a couple of co-workers. These deletions did not affect subgenomic RNAs (sgRNA) production downstream of ORF7a. Still, fragments associated with sgRNA of genes upstream of ORF7a showed a decrease in size when corresponding to samples with deletions. In silico analysis suggests that the deletions impair protein proper function; however, isolated viruses with partial deletion of ORF7a can replicate in culture cells similarly to wild-type viruses at 24 hpi, but with less infectious particles after 48 hpi. These findings on deleted ORF7a accessory protein gene, contribute to understanding SARS-CoV-2 phenotypes such as replication, immune evasion and evolutionary fitness as well insights into the role of SARS-CoV-2_ORF7a in the mechanism of virus-host interactions.

Insertion and deletion mutations preserved in SARS-CoV-2 variants.

The insertion/deletion (indel) mutation profiles of SARS-CoV-2 variants, including Omicron, remain unclear. We compared whole-genome sequences from various lineages and used preserved indels to infer the ancestral relationships between different lineages. Thirteen indel patterns from twelve sites were seen in ≥ 2 sequences; six of these sites were located in the N-terminal domain of the viral spike gene. Preserved indels in the coding regions were also identified in the non-structural protein 3 (Nsp3), Nsp6, and nucleocapsid genes. Seven of the thirteen indel patterns were specific to the Omicron variants, four of which were observed in BA.1, making it the most mutated variant. Other preserved indels observed in the Omicron variants were also seen in Alpha and/or Gamma, but not Delta, suggesting that Omicron is phylogenetically more proximal to Alpha. We demonstrated distinct profiles of preserved indels among SARS-CoV-2 variants and sublineages, suggesting the importance of indels in viral evolution.

Deletion of the s2m RNA Structure in the Avian Coronavirus Infectious Bronchitis Virus and Human Astrovirus Results in Sequence Insertions.

Coronaviruses infect a wide variety of host species, resulting in a range of diseases in both humans and animals. The coronavirus genome consists of a large positive-sense single-stranded molecule of RNA containing many RNA structures. One structure, denoted s2m and consisting of 41 nucleotides, is located within the 3' untranslated region (3' UTR) and is shared between some coronavirus species, including infectious bronchitis virus (IBV), severe acute respiratory syndrome coronavirus (SARS-CoV), and SARS-CoV-2, as well as other pathogens, including human astrovirus. Using a reverse genetic system to generate recombinant viruses, we investigated the requirement of the s2m structure in the replication of IBV, a globally distributed economically important Gammacoronavirus that infects poultry causing respiratory disease. Deletion of three nucleotides predicted to destabilize the canonical structure of the s2m or the deletion of the nucleotides corresponding to s2m impacted viral replication in vitro. In vitro passaging of the recombinant IBV with the s2m sequence deleted resulted in a 36-nucleotide insertion in place of the deletion, which was identified to be composed of a duplication of flanking sequences. A similar result was observed following serial passage of human astrovirus with a deleted s2m sequence. RNA modeling indicated that deletion of the nucleotides corresponding to the s2m impacted other RNA structures present in the IBV 3' UTR. Our results indicated for both IBV and human astrovirus a preference for nucleotide occupation in the genome location corresponding to the s2m, which is independent of the specific s2m sequence. IMPORTANCE Coronaviruses infect many species, including humans and animals, with substantial effects on livestock, particularly with respect to poultry. The coronavirus RNA genome consists of structural elements involved in viral replication whose roles are poorly understood. We investigated the requirement of the RNA structural element s2m in the replication of the Gammacoronavirus infectious bronchitis virus, an economically important viral pathogen of poultry. Using reverse genetics to generate recombinant IBVs with either a disrupted or deleted s2m, we showed that the s2m is not required for viral replication in cell culture; however, replication is decreased in tracheal tissue, suggesting a role for the s2m in the natural host. Passaging of these viruses as well as human astrovirus lacking the s2m sequence demonstrated a preference for nucleotide occupation, independent of the s2m sequence. RNA modeling suggested deletion of the s2m may negatively impact other essential RNA structures.

Expanding repertoire of SARS-CoV-2 deletion mutations contributes to evolution of highly transmissible variants.

The emergence of highly transmissible SARS-CoV-2 variants and vaccine breakthrough infections globally mandated the characterization of the immuno-evasive features of SARS-CoV-2. Here, we systematically analyzed 2.13 million SARS-CoV-2 genomes from 188 countries/territories (up to June 2021) and performed whole-genome viral sequencing from 102 COVID-19 patients, including 43 vaccine breakthrough infections. We identified 92 Spike protein mutations that increased in prevalence during at least one surge in SARS-CoV-2 test positivity in any country over a 3-month window. Deletions in the Spike protein N-terminal domain were highly enriched for these 'surge-associated mutations' (Odds Ratio = 14.19, 95% CI 6.15-32.75, p value = 3.41 × 10-10). Based on a longitudinal analysis of mutational prevalence globally, we found an expanding repertoire of Spike protein deletions proximal to an antigenic supersite in the N-terminal domain that may be one of the key contributors to the evolution of highly transmissible variants. Finally, we generated clinically annotated SARS-CoV-2 whole genome sequences from 102 patients and identified 107 unique mutations, including 78 substitutions and 29 deletions. In five patients, we identified distinct deletions between residues 85-90, which reside within a linear B cell epitope. Deletions in this region arose contemporaneously on a diverse background of variants across the globe since December 2020. Overall, our findings based on genomic-epidemiology and clinical surveillance suggest that the genomic deletion of dispensable antigenic regions in SARS-CoV-2 may contribute to the evasion of immune responses and the evolution of highly transmissible variants.

Remdesivir does not affect mitochondrial DNA copy number or deletion mutation frequency in aged male rats: A short report.

Remdesivir is a leading therapy in patients with moderate to severe coronavirus 2 (SARS-CoV-2) infection; the majority of whom are older individuals. Remdesivir is a nucleoside analog that incorporates into nascent viral RNA, inhibiting RNA-directed RNA polymerases, including that of SARS-CoV-2. Less is known about remdesivir's effects on mitochondria, particularly in older adults where mitochondria are known to be dysfunctional. Furthermore, its effect on age-induced mitochondrial mutations and copy number has not been previously studied. We hypothesized that remdesivir adversely affects mtDNA copy number and deletion mutation frequency in aged rodents. To test this hypothesis, 30-month-old male F333BNF1 rats were treated with remdesivir for three months. To determine if remdesivir adversely affects mtDNA, we measured copy number and mtDNA deletion frequency in rat hearts, kidneys, and skeletal muscles using digital PCR. We found no effects from three months of remdesivir treatment on mtDNA copy number or deletion mutation frequency in 33-month-old rats. These data support the notion that remdesivir does not compromise mtDNA quality or quantity at old age in mammals. Future work should focus on examining additional tissues such as brain and liver, and extend testing to human clinical samples.

Statistical analysis supports UTR (untranslated region) deletion theory in SARS-CoV-2.

It was noticed that the mortality rate of SARS-CoV-2 infection experienced a significant declination in the early stage of the epidemic. We suspect that the sharp deterioration of virus toxicity is related to the deletion of the untranslated region (UTR) of the virus genome. It was found that the genome length of SARS-CoV-2 engaged a significant truncation due to UTR deletion after a mega-sequence analysis. Sequence similarity analysis further indicated that short UTR strains originated from its long UTR ancestors after an irreversible deletion. A good correlation was discovered between genome length and mortality, which demonstrated that the deletion of the virus UTR significantly affected the toxicity of the virus. This correlation was further confirmed in a significance analysis of the genetic influence on the clinical outcomes. The viral genome length of hospitalized patients was significantly more extensive than that of asymptomatic patients. In contrast, the viral genome length of asymptomatic was considerably longer than that of ordinary patients with symptoms. A genome-level mutation scanning was performed to systematically evaluate the influence of mutations at each position on virulence. The results indicated that UTR deletion was the primary driving force in alternating virus virulence in the early evolution. In the end, we proposed a mathematical model to explain why this UTR deletion was not continuous.

Insertion-and-Deletion Mutations between the Genomes of SARS-CoV, SARS-CoV-2, and Bat Coronavirus RaTG13.

The evolutional process of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) development remains inconclusive. This study compared the genome sequences of severe acute respiratory syndrome coronavirus (SARS-CoV), bat coronavirus RaTG13, and SARS-CoV-2. In total, the genomes of SARS-CoV-2 and RaTG13 were 77.9% and 77.7% identical to the genome of SARS-CoV, respectively. A total of 3.6% (1,068 bases) of the SARS-CoV-2 genome was derived from insertion and/or deletion (indel) mutations, and 18.6% (5,548 bases) was from point mutations from the genome of SARS-CoV. At least 35 indel sites were confirmed in the genome of SARS-CoV-2, in which 17 were with ≥10 consecutive bases long. Ten of these relatively long indels were located in the spike (S) gene, five in nonstructural protein 3 (Nsp3) gene of open reading frame (ORF) 1a, and one in ORF8 and noncoding region. Seventeen (48.6%) of the 35 indels were based on insertion-and-deletion mutations with exchanged gene sequences of 7-325 consecutive bases. Almost the complete ORF8 gene was replaced by a single 325 consecutive base-long indel. The distribution of these indels was roughly in accordance with the distribution of the rate of point mutation rate around the indels. The genome sequence of SARS-CoV-2 was 96.0% identical to that of RaTG13. There was no long insertion-and-deletion mutation between the genomes of RaTG13 and SARS-CoV-2. The findings of the uneven distribution of multiple indels and the presence of multiple long insertion-and-deletion mutations with exchanged consecutive base sequences in the viral genome may provide insights into SARS-CoV-2 development. IMPORTANCE The developmental mechanism of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) remains inconclusive. This study compared the base sequence one-by-one between severe acute respiratory syndrome coronavirus (SARS-CoV) or bat coronavirus RaTG13 and SARS-CoV-2. The genomes of SARS-CoV-2 and RaTG13 were 77.9% and 77.7% identical to the genome of SARS-CoV, respectively. Seventeen of the 35 sites with insertion and/or deletion mutations between SARS-CoV-2 and SARS-CoV were based on insertion-and-deletion mutations with the replacement of 7-325 consecutive bases. Most of these long insertion-and-deletion sites were concentrated in the nonstructural protein 3 (Nsp3) gene of open reading frame (ORF) 1a, S1 domain of the spike protein, and ORF8 genes. Such long insertion-and-deletion mutations were not observed between the genomes of RaTG13 and SARS-CoV-2. The presence of multiple long insertion-and-deletion mutations in the genome of SARS-CoV-2 and their uneven distributions may provide further insights into the development of the virus.

Transmission cluster of COVID-19 cases from Uruguay: emergence and spreading of a novel SARS-CoV-2 ORF6 deletion.

BACKGROUND: Evolutionary changes in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) include indels in non-structural, structural, and accessory open reading frames (ORFs) or genes. OBJECTIVES: We track indels in accessory ORFs to infer evolutionary gene patterns and epidemiological links between outbreaks. METHODS: Genomes from Coronavirus disease 2019 (COVID-19) case-patients were Illumina sequenced using ARTIC_V3. The assembled genomes were analysed to detect substitutions and indels. FINDINGS: We reported the emergence and spread of a unique 4-nucleotide deletion in the accessory ORF6, an interesting gene with immune modulation activity. The deletion in ORF6 removes one repeat unit of a two 4-nucleotide repeat, which shows that directly repeated sequences in the SARS-CoV-2 genome are associated with indels, even outside the context of extended repeat regions. The 4-nucleotide deletion produces a frameshifting change that results in a protein with two inserted amino acids, increasing the coding information of this accessory ORF. Epidemiological and genomic data indicate that the deletion variant has a single common ancestor and was initially detected in a health care outbreak and later in other COVID-19 cases, establishing a transmission cluster in the Uruguayan population. MAIN CONCLUSIONS: Our findings provide evidence for the origin and spread of deletion variants and emphasise indels' importance in epidemiological studies, including differentiating consecutive outbreaks occurring in the same health facility.

From delta to Omicron: S1-RBD/S2 mutation/deletion equilibrium in SARS-CoV-2 defined variants.

Coronavirus-related Severe Acute Respiratory Syndrome (SARS-CoV) in 2002/2003, Middle-East Respiratory Syndrome (MERS-CoV) in 2012/2013, and especially the current 2019/2021 Severe Acute Respiratory Syndrome-2 (SARS-CoV-2) affected negatively the national health systems' endurance worldwide. SARS-Cov-2 virus belongs to lineage b of beta-CoVs demonstrating a strong phylogenetic similarity with BatCoVRaTG13 type. Spike (S) glycoprotein projections -consisting of two subunits S1/S2- provide a unique crown-like formation (corona) on virion's surface. Concerning their functional role, S1 represents the main receptor-binding domain (RBD), whereas S2 is involved in the virus-cell membrane fusion mechanism. On Nov 26th 2021, WHO designated the new SARS-CoV-2 strain - named Omicron, from letter ''όµÎ¹κρον'' in the Greek alphabet - as a variant of concern (B.1.1529 variant). Potentially this new variant is associated with high transmissibility leading to elevated infectivity and probably increased re-infection rates. Its impact on morbidity/mortality remains under investigation. In the current paper, analyzing and comparing the alterations of SARS-CoV-2 S RNA sequences in the defined variants (Alpha to Omicron), we observed some interesting findings regarding the S1-RBD/S2 mutation/deletion equilibrium that maybe affect and modify its activity.

A deletion in SARS-CoV-2 ORF7 identified in COVID-19 outbreak in Uruguay.

The analysis of genetic diversity in SARS-CoV-2 is the focus of several studies, providing insights into how the virus emerged and evolves. Most common changes in SARS-CoV-2 are single or point nucleotide substitutions; meanwhile, insertions and deletions (indels) have been identified as a less frequent source of viral genetic variability. Here, we report the emergence of a 12-nucleotide deletion in ORF7a, resulting in a 4-amino acid in-frame deletion. The Δ12 variant was identified in viruses from patients of a single outbreak and represents the first report of this deletion in South American isolates. Phylogenetic analysis revealed that Δ12 strains belong to the lineage B.1.1 and clustered separated from the remaining Uruguayan strains. The ∆12 variant was detected in 14 patients of this outbreak by NGS sequencing and/or two rapid and economic methodologies: Sanger amplicon sequencing and capillary electrophoresis. The presence of strong molecular markers as the deletion described here are useful for tracking outbreaks and reveal a significant aspect of the SARS-CoV-2 evolution on the robustness of the virus to keep its functionality regardless loss of genetic material.

Structural analysis of SARS-Cov-2 nonstructural protein 1 polymorphisms found in the Brazilian Amazon.

The coronavirus disease COVID-19 has been the cause of millions of deaths worldwide. Among the SARS-CoV-2 proteins, the non-structural protein 1 (NSP1) has great importance during the virus infection process and is present in both alpha and beta-CoVs. Therefore, monitoring of NSP1 polymorphisms is crucial in order to understand their role during infection and virus-induced pathogenicity. Herein, we analyzed how mutations detected in the circulating SARS-CoV-2 in the population of the city of Manaus, Amazonas state, Brazil could modify the tertiary structure of the NSP1 protein. Three mutations were detected in the SARS-CoV-2 NSP1 gene: deletion of the amino acids KSF from positions 141 to 143 (delKSF), SARS-CoV-2, lineage B.1.195; and two substitutions, R29H and R43C, SARS-CoV-2 lineage B.1.1.28 and B.1.1.33, respectively. The delKSF was found in 47 samples, whereas R29H and R43C were found in two samples, one for each mutation. The NSP1 structures carrying the mutations R43C and R29H on the N-terminal portion (e.g. residues 10 to 127) showed minor backbone divergence compared to the Wuhan model. However, the NSP1 C-terminal region (residues 145 to 180) was severely affected in the delKSF and R29H mutants. The intermediate variable region (residues 144 to 148) leads to changes in the C-terminal region, particularly in the delKSF structure. New investigations must be carried out to analyze how these changes affect NSP1 activity during the infection. Our results reinforce the need for continuous genomic surveillance of SARS-CoV-2 to better understand virus evolution and assess the potential impact of the viral mutations on the approved vaccines and future therapies.

Multiple Occurrences of a 168-Nucleotide Deletion in SARS-CoV-2 ORF8, Unnoticed by Standard Amplicon Sequencing and Variant Calling Pipelines.

Genomic surveillance of the SARS-CoV-2 pandemic is crucial and mainly achieved by amplicon sequencing protocols. Overlapping tiled-amplicons are generated to establish contiguous SARS-CoV-2 genome sequences, which enable the precise resolution of infection chains and outbreaks. We investigated a SARS-CoV-2 outbreak in a local hospital and used nanopore sequencing with a modified ARTIC protocol employing 1200 bp long amplicons. We detected a long deletion of 168 nucleotides in the ORF8 gene in 76 samples from the hospital outbreak. This deletion is difficult to identify with the classical amplicon sequencing procedures since it removes two amplicon primer-binding sites. We analyzed public SARS-CoV-2 sequences and sequencing read data from ENA and identified the same deletion in over 100 genomes belonging to different lineages of SARS-CoV-2, pointing to a mutation hotspot or to positive selection. In almost all cases, the deletion was not represented in the virus genome sequence after consensus building. Additionally, further database searches point to other deletions in the ORF8 coding region that have never been reported by the standard data analysis pipelines. These findings and the fact that ORF8 is especially prone to deletions, make a clear case for the urgent necessity of public availability of the raw data for this and other large deletions that might change the physiology of the virus towards endemism.

A SARS-CoV-2 mutant from B.1.258 lineage with ∆H69/∆V70 deletion in the Spike protein circulating in Central Europe in the fall 2020.

SARS-CoV-2 mutants carrying the ∆H69/∆V70 deletion in the amino-terminal domain of the Spike protein emerged independently in at least six lineages of the virus (namely, B.1.1.7, B.1.1.298, B.1.160, B.1.177, B.1.258, B.1.375). We analyzed SARS-CoV-2 samples collected from various regions of Slovakia between November and December 2020 that were presumed to contain B.1.1.7 variant due to drop-out of the Spike gene target in an RT-qPCR test caused by this deletion. Sequencing of these samples revealed that although in some cases the samples were indeed confirmed as B.1.1.7, a substantial fraction of samples contained another ∆H69/∆V70 carrying mutant belonging to the lineage B.1.258, which has been circulating in Central Europe since August 2020, long before the import of B.1.1.7. Phylogenetic analysis shows that the early sublineage of B.1.258 acquired the N439K substitution in the receptor-binding domain (RBD) of the Spike protein and, later on, also the deletion ∆H69/∆V70 in the Spike N-terminal domain (NTD). This variant was particularly common in several European countries including the Czech Republic and Slovakia but has been quickly replaced by B.1.1.7 early in 2021.

Antibody evasion by the P.1 strain of SARS-CoV-2.

Terminating the SARS-CoV-2 pandemic relies upon pan-global vaccination. Current vaccines elicit neutralizing antibody responses to the virus spike derived from early isolates. However, new strains have emerged with multiple mutations, including P.1 from Brazil, B.1.351 from South Africa, and B.1.1.7 from the UK (12, 10, and 9 changes in the spike, respectively). All have mutations in the ACE2 binding site, with P.1 and B.1.351 having a virtually identical triplet (E484K, K417N/T, and N501Y), which we show confer similar increased affinity for ACE2. We show that, surprisingly, P.1 is significantly less resistant to naturally acquired or vaccine-induced antibody responses than B.1.351, suggesting that changes outside the receptor-binding domain (RBD) impact neutralization. Monoclonal antibody (mAb) 222 neutralizes all three variants despite interacting with two of the ACE2-binding site mutations. We explain this through structural analysis and use the 222 light chain to largely restore neutralization potency to a major class of public antibodies.

Angiotensin-converting enzyme-1 gene insertion/deletion polymorphism may be associated with COVID-19 clinical severity: a prospective cohort study.

BACKGROUND: Angiotensin-converting enzyme (ACE) insertion/deletion (I/D) polymorphism may play a role in the pathogenesis of coronavirus-19 disease (COVID-19). OBJECTIVES: Investigate the relationship between ACE I/D polymorphism and the clinical severity of COVID-19. DESIGN: Prospective cohort study. SETTING: Tertiary care hospital. PATIENTS AND METHODS: The study included COVID-19 patients with asymptomatic, mild, and severe disease with clinical data and whole blood samples collected from 1 April 2020 to 1 July 2020. ACE I/D genotypes were determined by polymerase chain reaction and agarose gel electrophoresis. MAIN OUTCOME MEASURE: ACE DD, DI and II genotypes frequencies. SAMPLE SIZE: 90 cases, 30 in each disease severity group. RESULTS: Age and the frequency of general comorbidity increased significantly from the asymptomatic disease group to the severe disease group. Advanced age, diabetes mellitus and presence of ischemic heart disease were independent risk factors for severe COVID-19 [OR and 95 % CI: 1.052 (1.021-1.083), 5.204 (1.006-26.892) and 5.922 (1.109-31.633), respectively]. The ACE II genotype was the dominant genotype (50%) in asymptomatic patients, while the DD genotype was the dominant genotype (63.3 %) in severe disease. The ACE II geno-type was protective against severe COVID-19 [OR and 95% CI: .323 (.112-.929)]. All nine patients (8.9%) who died had severe disease. CONCLUSIONS: The clinical severity of COVID-19 infection may be associated with the ACE I/D polymorphism. LIMITATIONS: Small sample size and single center. CONFLICT OF INTEREST: None.

SARS-CoV-2's claimed natural origin is undermined by issues with genome sequences of its relative strains: Coronavirus sequences RaTG13, MP789 and RmYN02 raise multiple questions to be critically addressed by the scientific community.

RaTG13, MP789, and RmYN02 are the strains closest to SARS-CoV-2, and their existence came to light only after the start of the pandemic. Their genomes have been used to support a natural origin of SARS-CoV-2 but after a close examination all of them exhibit several issues. We specifically address the presence in RmYN02 and closely related RacCSxxx strains of a claimed natural PAA/PVA amino acid insertion at the S1/S2 junction of their spike protein at the same position where the PRRA insertion in SARS-CoV-2 has created a polybasic furin cleavage site. We show that RmYN02/RacCSxxx instead of the claimed insertion carry a 6-nucleotide deletion in the region and that the 12-nucleotide insertion in SARS-CoV-2 remains unique among Sarbecoviruses. Also, our analysis of RaTG13 and RmYN02's metagenomic datasets found unexpected reads which could indicate possible contamination. Because of their importance to inferring SARS-CoV-2's origin, we call for a careful reevaluation of RaTG13, MP789 and RmYN02 sequencing records and assembly methods.

The SARS-CoV-2 Transcriptome and the Dynamics of the S Gene Furin Cleavage Site in Primary Human Airway Epithelia.

The spike (S) polypeptide of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) consists of the S1 and S2 subunits and is processed by cellular proteases at the S1/S2 boundary that contains a furin cleavage site (FCS), 682RRAR↓S686 Various deletions surrounding the FCS have been identified in patients. When SARS-CoV-2 propagated in Vero cells, it acquired deletions surrounding the FCS. We studied the viral transcriptome in Vero cell-derived SARS-CoV-2-infected primary human airway epithelia (HAE) cultured at an air-liquid interface (ALI) with an emphasis on the viral genome stability of the FCS. While we found overall the viral transcriptome is similar to that generated from infected Vero cells, we identified a high percentage of mutated viral genome and transcripts in HAE-ALI. Two highly frequent deletions were found at the FCS region: a 12 amino acid deletion (678TNSPRRAR↓SVAS689) that contains the underlined FCS and a 5 amino acid deletion (675QTQTN679) that is two amino acids upstream of the FCS. Further studies on the dynamics of the FCS deletions in apically released virions from 11 infected HAE-ALI cultures of both healthy and lung disease donors revealed that the selective pressure for the FCS maintains the FCS stably in 9 HAE-ALI cultures but with 2 exceptions, in which the FCS deletions are retained at a high rate of >40% after infection of ≥13 days. Our study presents evidence for the role of unique properties of human airway epithelia in the dynamics of the FCS region during infection of human airways, which is likely donor dependent.IMPORTANCE Polarized human airway epithelia at an air-liquid interface (HAE-ALI) are an in vitro model that supports efficient infection of SARS-CoV-2. The spike (S) protein of SARS-CoV-2 contains a furin cleavage site (FCS) at the boundary of the S1 and S2 domains which distinguishes it from SARS-CoV. However, FCS deletion mutants have been identified in patients and in vitro cell cultures, and how the airway epithelial cells maintain the unique FCS remains unknown. We found that HAE-ALI cultures were capable of suppressing two prevalent FCS deletion mutants (Δ678TNSPRRAR↓SVAS689 and Δ675QTQTN679) that were selected during propagation in Vero cells. While such suppression was observed in 9 out of 11 of the tested HAE-ALI cultures derived from independent donors, 2 exceptions that retained a high rate of FCS deletions were also found. Our results present evidence of the donor-dependent properties of human airway epithelia in the evolution of the FCS during infection.

Genome sequencing of SARS-CoV-2 in a cohort of Egyptian patients revealed mutation hotspots that are related to clinical outcomes.

BACKGROUND: Severe acute respiratory syndrome-2 (SARS-CoV-2) exhibits a broad spectrum of clinical manifestations. Despite the fact that SARS-CoV-2 has slower evolutionary rate than other coronaviruses, different mutational hotspots have been identified along the SARS-CoV-2 genome. METHODS: We performed whole-genome high throughput sequencing on isolates from 50 Egyptian patients to see if the variation in clinical symptoms was related to mutations in the SARS-CoV-2 genome. Then, we investigated the relationship between the observed mutations and the clinical characteristics of the patients. RESULTS: Among the 36 most common mutations, we found two frameshift deletions linked to an increased risk of shortness of breath, a V6 deletion in the spike glycoprotein's signal peptide region linked to an increased risk of fever, longer fever duration and nasal congestion, and L3606-nsp6 deletion linked to a higher prevalence of cough and conjunctival congestion. S5398L nsp13-helicase was linked to an increased risk of fever duration and progression. The most common mutations (241, 3037, 14,408, and 23,403) were not linked to clinical variability. However, the E3909G-nsp7 variant was more common in children (2-13 years old) and was associated with a shorter duration of symptoms. The duration of fever was significantly reduced with E1363D-nsp3 and E3073A-nsp4. CONCLUSIONS: The most common mutations, D614G/spike-glycoprotein and P4715L/RNA-dependent-RNA-polymerase, were linked to transmissibility regardless of symptom variability. E3909G-nsp7 could explain why children recover so quickly. Nsp6-L3606fs, spike-glycoprotein-V6fs, and nsp13-S5398L variants may be linked to clinical symptom worsening. These variations related to host-virus interactions might open new therapeutic avenues for symptom relief and disease containment.